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LAPLACE OPERATORS WITH EIGENFUNCTIONS 
WHOSE NODAL SET IS A KNOT 

ALBERTO ENCISO, DAVID HARTLEY, AND DANIEL PERALTA-SALAS 


Abstract. We prove that, given any knot 7 in a compact 3-manifold M, 
there exists a Riemannian metric on M such that there is a complex-valued 
eigenfunction u of the Laplacian, corresponding to the first nontrivial eigen¬ 
value, whose nodal set has a connected component given by 7 . Higher 

dimensional analogs of this result will also be considered. 


1. Introduction 

Let M be a closed manifold of dimension 3 endowed with a smooth Riemannian 
metric g. The eigenfunctions of M satisfy the equation 

Auk = —XkUk , 

where 0 = Aq < Ai ^ A2 ^ . are the eigenvalues of M and A is the Laplace 

operator of the manifold. As is well known, the zero set u^^(O) is called the nodal 
set of the eigenfunction. 

Our goal in this paper is to establish the existence of knotted nodal sets in eigen¬ 
functions of the Laplacian on a Riemannian 3-manifold. Specifically, the starting 
point of the present paper is the following 

Question. Let 7 be a (possibly knotted) curve in M. Is there a complex-valued 
eigenfunction u of the Laplacian whose nodal set is diffeomorphic to 7? 

Before getting into the background of this question, it is worth discussing why 
the eigenfunction u is assumed complex-valued. The idea is that, by the maximum 
principle, the nodal set of a real-valued eigenfunction (i.e., u)r^(O)) cannot have a 
connected component of codimension greater than or equal to 2, so there cannot be 
a real eigenfunction of the Laplacian whose nodal set has a component diffeomorphic 
to the curve 7. A way to bypass this obstruction is to take a complex function of 
the form, say, u := ui + iu2, which will indeed be an eigenfunction provided that 
Ai = A2. Notice that u“^(0) = u)”^(0) fl u^^(O). Hence the above question is in 
fact a problem about Laplacians with degenerate eigenvalues. 

The context of this question is the study of the nodal sets of the eigenfunctions 
of the Laplacian in a compact Riemannian manifold, which is a classical topic in 
geometric analysis with a number of important open problems [201 HI]. In a way, 
this question addresses the flexibility of low-energy Laplace eigenfunctions, as we 
shall explain next. 

It is well known that high-energy eigenvalues and eigenfunctions of the Laplacian 
are rather rigid, as illustrated by Weyl’s asymptotic formula for the eigenvalues and 
by Yau’s conjecture concerning the Hausdorff measure of the eigenfunction’s nodal 
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sets as the eigenvalue tends to infinity, which was proved for real-analytic metrics 
in [6]. Indeed, the measure-theoretic properties of the nodal set of high-energy 
eigenfunctions have been a hot topic of research for decades umin] ; a recent 
account of the subject can be found in m- 

On the contrary, the situation for low-energy eigenvalues and eigenfunctions 
is much more flexible. The eigenvalue case was considered by Colin de Verdiere 
in [^, where it was shown that there are metrics whose first N eigenvalues can be 
prescribed a priori. Concerning nodal sets, the problem of constructing a metric 
on M for which the nodal set of the first eigenfunction has a prescribed 

connected component was solved in m (see also [IS US] for the two dimensional 
case). 

The question that we have started with addresses the possibility of prescribing 
a component of the intersection of nodal sets 

M-i(O) = uj“^(0) nM^^(O), 

subject to the additional constraint that the eigenvalues Ai and A2 should coincide. 
This constraint is important, since it is well known that for a generic metric the 
eigenvalues are all simple | 19 ] . In a way, this problem is the metric counterpart of 
Berry’s conjecture [ 3 ], recently proved in [ 7 ], which states that there should be a 
smooth (real) potential for which the Schrddinger operator — A -|- 1 / in admits 
a complex-valued eigenfunction whose nodal set has a connected component of any 
given knot type. 

Our main result in this paper provides an affirmative answer to the above ques¬ 
tion. Specifically, we show that given a smooth closed curve 7 (without self¬ 
intersections, but possibly knotted), there is a metric for which u“^(0) has a 
connected component given by 7, with u := ui + m2 corresponding to the first 
eigenvalue Ai = A2, which has multiplicity 2 . Moreover, the knot 7 is structurally 
stable in the sense that any small enough perturbation of u (in the norm with 
k ^ \) has a connected component in its nodal set that is a small perturbation 
of 7 (i.e., the image of 7 under a smooth diffeomorphism of M that is close to the 
identity in the norm). We are thus led to the following: 

Theorem 1.1. Let 7 6e o (possibly knotted) closed curve in a closed 3-manifold 
M. Then there exists a Riemannian metric g on M such that its first nontrivial 
eigenvalue has multiplicity 2 (i.e., Ai = A 2 ^ Asj and j is a connected component 
of the nodal set u“^(0) of the complex-valued eigenfunction u = ui-\-iu 2 . Moreover, 
') is a structurally stable component of the nodal set. 

Equivalently, the theorem asserts that the closed curve 7 is a connected compo¬ 
nent of the intersection of the nodal sets of the real-valued eigenfunctions uj"^(0) fl 
M^^(O) for any basis of the eigenspace corresponding to the first nontrivial eigen¬ 
value. 

The proof of this result starts off with the construction of a metric g that is of 
order e in the complement of a tubular neighborhood fl of the curve 7. If e is small 
enough, the eigenfunctions ui and U 2 are close in the set fl to the first nontrivial 
Neumann eigenfunctions vi and V 2 of this set, so in principle the nodal set of the 
former can be controlled in terms of the latter using Thom’s isotopy theorem. To 
develop this strategy there are two difficulties one has to deal with. First, the 
intersection uj~^( 0 ) fl u^^(O) must be transverse in order to apply Thom’s isotopy 
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theorem. Second, the eigenvalue Ai must be degenerate, which is a non-generic 
property according to Uhlenbeck’s theorem m- These problems are overcome by 
constructing very carefully the metric g in successive steps. 

Since the nodal set of the complex-valued eigenfunction u = ui + iu 2 is the 
intersection of the nodal sets u^^( 0 )n'u^^( 0 ), the question that we have started our 
paper with is simply the easiest case of whether we can realize a set of codimension 
m ^ 2 in a compact d-manifold as the intersection of the nodal sets 

u^i(o)n---nu-^(o) 

of m real-valued eigenfunctions having the same energy Ai = • • • = Am (of mul¬ 
tiplicity m). The full problem will be tackled in Section O where we provide a 
higher-dimensional version of Theorem 11.11 iTheorem 15.11) . This is a significant 
generalization of cni, but its statement is a little more technical than that of The¬ 
orem [IT] 

The content of this paper is as follows. The proof of Theorem 11.11 is presented 
in Section [21 although the proofs of two technical results. Propositions 12.II and 12.21 
are relegated to Sections [3] and lU respectively. Finally, in Section [5] we state and 
prove a higher dimensional analog of Theorem 11.11 


2. Proof of the main theorem 

The proof of Theorem 11.11 is divided in four steps. In Step 1 we introduce an 
appropriate metric go whose explicit expression in a tubular neighborhood n of the 
curve 7 allows us to compute, in closed form, the first Neumann eigenvalues and 
eigenfunctions of the domain fl. In Step 2 we rescale the metric go with a smooth 
factor that is of order e in the complement of fl, and show that the eigenvalues 
and eigenfunctions of M with this metric tend to the Neumann eigenvalues and 
eigenfunctions of 11 as e \ 0 (Proposition 12.II) . In Step 3 we prove that the metric 
introduced in the previous step can be deformed so that its first nontrivial eigenvalue 
has multiplicity 2 and its corresponding eigenspace is close to an eigenspace of 
Neumann eigenfunctions of H (Proposition 12.21) . Finally, in Step 4 we use Thom’s 
isotopy theorem to control the intersection of the nodal sets of the eigenfunctions 
obtained before. 


Step 1: Neumann eigenfunctions of a certain bounded domain. Let us consider 
a tubular neighborhood 12 of the curve 7 , which is diffeomorphic to the product 
X because the normal bundle of a knot is always trivial [15]. Hence we will 
parametrize the set H using coordinates z = {zi,Z 2 ) G and s G M/Z. Here 
denotes the unit two-dimensional disk. 


We shall take a metric go on M such that its restriction to H can be written in 
these coordinates as 


5 o|n := A ^ ds^ + |dzp = A ^ ds^ + dr'^ + dO'^ . 


Here (r, 9) are the polar coordinates in the disk associated to z and T is a large 
enough constant. Specifically, we assume that 


( 2 . 1 ) 


A > 


47r^ 


where is the first nontrivial Neumann eigenvalue of the unit 2-disk. Notice 
that the coordinate s essentially corresponds to a parametrization of the curve 7 
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proportional to arc-length and that the variables z describe a disk orthogonal to 
the curve. 

Let us consider the Neumann eigenvalue problem 

Avk = —fJ-kVk in ^ , d^Vk = 0 on dO ., 

where A stands for the Laplacian associated with the metric go- By separation of 
variables and using the condition (EH), we infer that the Neumann eigenvalues fik 
of the domain with respect to the metric go are then go = 0 , 

= ^2 = , 

and /ifc > g-D 2 for fc ^ 3. Hence the eigenspace of gi is two-dimensional, with a 
particular orthogonal basis being 

= Ji(v 7 ht) COS 0 , V 2 = Ji{y/g^r) sinO. 

Here Ji denotes the Bessel function. 

Since Ji{y/JIir) does not vanish for 0 < r < 1 and the eigenfunctions behave in 
a neighborhood of the origin as 

cose + 0{r^) = + OUz^ , = ^z^ + Odzp), 

it is elementary that, as the eigenfunctions are only defined in the set \z\ < 1 , the 
intersection 

is the set z = 0 , that is, the curve 7 (which is identified via the coordinates with 
the set §1 X {0}). Moreover, the behavior of the Neumann eigenfunctions near the 
set z = 0 shows that their intersection is transverse in the sense that 

(2.2) rank(Vi;i(x), Vw 2 (ic)) = 2 

for all X G i;fi(0) n p^ 1(0). These properties obviously hold for any other basis 
{vi,V 2 } of the eigenspace corresponding to gi. 


Step 2: Convergence of eigenvalues and eigenfunctions. Let us take a one-parameter 
family of smooth, positive, uniformly bounded functions ft on M such that ft{x) = e 
in the complement := M\H and ft{x) = 1 if a: G H and the distance between x 
and dfl, as measured with the metric go, is larger than 1/t. Here both t and e are 
positive real numbers, and ft is assumed to be analytic in t and e. Notice that the 
smooth functions ft converge pointwise to the discontinuous function 


fix) 


e if a; G , 
1 if a; G H , 


as t —>■ 00 . 

Let us now define the one-parameter family of smooth metrics gt ■= ft go- In the 
following proposition we consider the behavior of the eigenvalues Xk^t and eigen¬ 
functions Ukg of the metric gt as the parameter e tends to zero and t tends to 
infinity. Specifically, we are interested in comparing these quantities with the Neu¬ 
mann eigenvalues gk and eigenfunctions Vk of the domain H, which were introduced 
in Step 1: 
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Proposition 2.1. There is a basis of the eigenfunctions Uk,t of the Laplacian on M 
with the metric gt and a basis Vk of the Neumann eigenfunctions of the domain 
such that 

Pk — c(l) , 

\\uk,t — Vk\\m{Q) = o(l). 

Here we are denoting by o(l) quantities that tend to zero (not uniformly in k) as 
e \ 0 and t ^ oo. 

The proof of this result is presented in Section |3l and uses ideas introduced by 
Colin de Verdiere [4]. 

Step 3: Collapse of eigenvalues. From Proposition [2T] we infer that the eigenvalues 
Ai^t and X 2 U are very close to the first Neumann eigenvalue provided that t 
is large enough and e is sufficiently small. Uhlenbeck’s theorem m implies that, 
generically, Ai^t ^ A 2 ,t, so generally if we want to get a sequence of eigenvalues of 
multiplicity 2 converging to we need to deform the family of metrics gt- In 
the following proposition we show that such a deformation exists. The statement 
of the proposition is in terms of an e-dependent smooth metric g^, e > 0, whose 
eigenvalues will be denoted by Xk,e- 

Proposition 2.2. There exists a positive e-dependent function h^ on M of class 
C°°, a basis of the eigenfunctions Uk,e of the Laplacian on M with the metric 
:= h^gr, where T := e~^, and a basis Vk of the Neumann eigenfunctions of the 
domain Tl (with the metric go) such that: 

(i) Ai^e has multiplicity 2 and Ai_e = A 2 ,e = 

(ii) limeN^o ||ui.£ - Pi||Hi(n) = lime\o \\u2,e - 11211 hi( n) = 0 . 

Moreover, limeN^o ll^e ~ 1 ||c‘(m) = 0 for any integer 1. 

The proof of this result is in Section 0] and exploits the fact that when e \ 0 the 
eigenvalues Ai,t and A 2 ,t are very close according to Proposition 12.11 (recall that 
gi = fj .2 = Md^)- Then a Hadamard-type formula for the variation of eigenvalues 
with respect to metric perturbations allows us to use an implicit function theorem 
argument to construct the function he- 

Step 4- Transverse intersection of nodal sets. Proposition 12.21 ensures that, for 
any fixed <5 > 0 , one can then choose a small e so that the eigenvalue Ai_e of the 
Laplacian computed with the metric has multiplicity 2, and that there is a basis 
{ui^e,U 2 ,e} of the corresponding eigenspace such that 

(2-3) ||Ml,e — IIl||Hi(n) < (5 , ||M 2 ,e — U 2 ||Hi(n) < ^ ■ 

Here {f 1 , 112 } is a basis of Neumann eigenfunctions of the domain LI. 

Since the conformal factor he defining the metric in Proposition 12.21 is arbi¬ 
trarily close to 1 as e 0 , the definition of the metric gt implies that 

Wle — 5o||c‘(S) < Ci 5 

for any integer I and any compact subset S G LI, where the constant Ci depends 
on I and S. Therefore, standard elliptic estimates allow us to promote the 
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bound (j2.3[l to a bound: 

(2-4) ||ui_£ — fill C1(5) < (7(5 , ||M 2 ,e — U 2 ||ci(S) < ^"(5 . 

For small enough 5 (and hence small e), the transversality property (12.21) allows 
us to apply Thom’s isotopy theorem [1] Section 20.2] to conclude that the nodal 
set u]((g( 0 ) n has a connected component in fl diffeomorphic to 

7 = ^;-i(0)nu2'\0), 

and that the corresponding diffeomorphism : M ^ M can be chosen arbitrarily 
close to the identity in C^{M) and different from the identity only in a small 
neighborhood of 7 . 

Finally, let us consider the pulled-back metric 

9 ■■= K9e 

and call 

Ui := o , U2 := U2,e => 

the associated eigenfunctions, whose corresponding eigenvalue Ai := Ai_e has mul¬ 
tiplicity 2. Accordingly, the nodal set w]"^(0) D u^^(O) has a connected component 
given by the curve 7 , which is structurally stable because Eq. (lO) and the 
estimates (12.41) imply 

rank(Vui(x), VM 2 (a:)) =2 
for all X G 7 . The theorem then follows. 

Remark 2.3. The transversality property (12.21) of the intersection of the nodal sets 
z;]"^(0) and u^^(O) is key in the proof, as in particular it allows us to use Thom’s 
isotopy theorem in Step 4. 


3. Proof of Proposition 12.11 


Since the family of smooth functions ft converges pointwise to the discontinuous 
function / as t —>■ 00 , it is convenient to work first with the discontinuous metric 


9e ■= fgo 


ego if X G , 

go if X G n , 


and later we will prove the desired result for the metric gt provided that t is large 
enough. 

To define the spectrum of the Laplacian associated to the metric we use the 
quadratic form 

Qe{(p) ■■= [ \d(p\ldVe= [ Idifl'^ + e^ [ |d(pp, 

JM Jn 

and the natural norm corresponding to the metric g^- 
||(^||^= / <p^dV,= [ (^^ + ei / 

JM JQ. Jvt‘ 

Here the subscripts e refer to quantities computed with respect to the metric and 
we are omitting the subscripts (and indeed the measure in the integrals) when the 
quantities correspond to the reference metric go- As is well known, the domain of 
the quadratic form (5e can be taken to be the Sobolev space (Recall that. 
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M being compact, this Sobolev space is independent of the smooth metric one uses 
to define it). 

By the min-max principle, the fc*** eigenvalue Xk,e of this quadratic form is 


(3.1) 




inf max 
w&Wk ¥>ew\{o} 




where Wk stands for the set of {k + l)-dimensional linear subspaces of H^{M) and 


9e(‘/5) := 




is the Rayleigh quotient. The fc**' eigenfunction Uk,€ is then a minimizer of the 
above variational problem for Afc,£, in the sense that any subspace that minimizes 
the variational problem can be written as span{uo,e,..., Uk,e}- 

In what follows we fix an integer feg- us first prove that the eigenvalues Xk^e 
are almost upper bounded by fj,k for k ^ kf) in the sense that 

(3.2) Xk,e ^ Mfc + 0(£^)) 


where 0(e™) stands for a quantity that is bounded by Ce™. For this purpose, 
we consider the harmonic extension to the whole manifold M of the Neumann 
eigenfunctions Vk of the domain fl: 

1^ Vk m il, 

where Vk is the solution to the boundary value problem 


ADfc = 0 in Dfclan = wfelan . 


Recall that A is the Laplacian computed with the reference metric go- Throughout 
this section we shall normalize the Neumann eigenfunctions Vk so that they define 
an orthonormal basis of L^{Vl). Standard elliptic estimates then give us the bound 

(3.3) \dvl\ + ^ ’ 

where C is a constant independent of e and k. 

For any linear combination (p,(p & span{'(/;o,..., V’fe}, with k ^ ko, using Eq. (13.31) 
we obtain that 


ipifi dVe = 


M 


+ ip(p 

Q Jq‘= 


= / (pif + 0{e^){gko+^)\\‘p\\L^n)\\<p\\L^Q), 


= (^(^ + 0(e2)||(/3||i2(n)||<p||L2(o), 

Jn 


where, to pass to the last line, we have absorbed the constant Uko +1 in the 0 (e 2 ) 
term. Accordingly, for sufficiently small e, the linear space 


span{i/io,.. ■,'ipk}, 
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is a (fc + l)-dimensional subspace of H^{M) provided that k ^ ko. In a similar 
manner we have 


Qe{^) = [ / \d(pf 

Jn 

~In ■ 


Therefore, for k ^ kg we can estimate Xk^e as 

Qeiv’) 


Afc,£ ^ max 

ipe{ipo,---,'ipk} 


^ max 
ipe{ipo,---,'<Pk} 


= max 
xe{^o,---,vk} 


Mhm 

ilh*!l + 0(e4) 

iLlli-, 1 ,, 


—Mfe + 0{e ^), 

which proves the upper bound (|3.2I) . 

Now we proceed to prove the lower bound 

(3.4) Afc_e ^ ^1 + 0(e2 fik 

for all k ^ kg. To obtain this bound we split {M) into two closed linear subspaces 
Ui := {ip G H^{M) : Av? = 0 in , 

H2 := [p G H\M) : p\nc G p\n = 0}, 

so that H^{M) = "Hi © 7^2 and 

[ 5o(VvJi,Vy)2) = 0, 

JM 

for all pj G TLj. We will write pj for the component of a function p G H^{M) in 

dir 

We hrst note that by choosing a sufficiently small e, the upper bound (13.21) 
implies that Xk,t < Pko + 1 for all fc ^ fcg, which allows us to reduce the min-max 
formula dSU for fc ^ fco to 

(3.5) Xk,e = inf max q^{p ), 

IVGW' ¥=elV\{0} 

where is the set of (fc + l)-dimensional subspaces of H^{M) such that 

q^{p) < pko+i 

for all nonzero p & W C. >V(,. 

A simple computation shows that 

(3.6) 11^211? = 0(e)||^||^ 

for any p & W C. with k Pi kQ. Indeed, this relation follows from the estimate 


IIV’ll? 11^211^ 


In- i’l 


e 11^11 = 
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where the e-independent constant An^ is the first Dirichlet eigenvalue of In 
particular, Eq. (j3.6l) implies that 

(3.7) 11^11? = (l + 0(e^))||<Pl||^ 

Proceeding as in the proof of Eq. (13.3|) and using the bound (13.7|) , we obtain that 
for G IE C with k < ko, 

f f \d^i\^^c [ 

Jn<= Jn<= Jn 


Combining these two results we conclude that, 


=(^1 + 0(6 = )) 


/m M'T’ile + /m dV^ 

WviWl 


^(l + 0 (e^)) 

= (l + 0(e^)) 

for any ip & W C with k ^ k^. 
have 


In\dVi\^ 

In In- 

fn^Pf 

Therefore by the min-max principle (13.51) we 


Afc,e ^ 




("l + 0 (e 2 )') inf max 
V / WGW' VJgW\{0} 


/nMy’iP 

fn^Pf 


+ 0(€^)'j Pk, 


for all k ^ ko, which establishes the lower bound (13.41) . 
The bounds (13.2p and (13.4p imply that 


(3.8) lim ||Afc,e - Mfcll = 0 

e\,0 

for all k ^ ko. Since the metric is equal to in the domain and the eigenvalues 
control the norm, it is standard that the convergence of eigenvalues (13.81) implies 
that for k ^ kf), there exists a basis of eigenfunctions Uk^e and a basis of Neumann 
eigenfunctions Vk such that 


(3.9) 


Im II- Pfe 11^1(0) =0. 


To finish the proof of the proposition, recall that the family of smooth metrics 
gt converges pointwise to g^. It is then well known (see e.g. [2]) that the eigenvalue 
Xk,t converges to Xk,e and that there is a basis of eigenfunctions Uk,t converging 
in H^{M) to a basis of eigenfunctions Ufc,£ as t —>■ oo. Combining this result with 
equations (|3.8[) and (13.91) . the proposition follows. 


4. Proof of Proposition 12.21 


In order to construct the conformal factor we introduce an auxiliary smooth 
function hg depending on two real parameters s = (si,S 2 ) G that has the form 

U _ „2sil/'l+2s2'i/’2 

ILc — C , 
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where '0i and ■02 are smooth functions with zero normal derivative on dfl that will 
be fixed later. We denote by Xk,e,s the eigenvalues of the smooth metric 

9e,s ■= hsgr , 

where T := and by Uk,e,s an orthonormal basis of the corresponding eigenfunc¬ 
tions. By construction we have that is analytic in the parameters e > 0 and 
s £ and that g^^o = gr- We also denote by gk,s the Neumann eigenvalues of 
the Laplacian associated with the smooth metric 

go,s ■= hsgo 

in the domain fl, and by Vk,B an orthonormal basis of the corresponding eigen¬ 
functions. By the Kato-Rellich theorem, since the Laplacian for the metric ge,s is 
analytic in the parameters s £ and e > 0, it is well known that we can take 
the families of eigenvalues Xk,e,s and eigenfunctions Uk,e,s to vary analytically in s 
and e (this result is independent of the multiplicities of the eigenvalues). The same 
theorem also implies that the Neumann eigenvalues fik,s and eigenfunctions Vk,s are 
analytic in s, and in particular 


(4.1) 


lim \ fik,s - gk\ = 0 , 

s —>-0 


lim - 'yfcllffi(n) = 0. 


To construct the desired function 
function theorem to the map A : R+ x R 


we will apply a version of the implicit 
' —7> R^ given by 


A(e,s) := 


(Al,e,s — , A2,e,s — ) if £ > 0 , 

(iti.s - MD2,/r2.s -/iD^) ife = 0. 


From the analytic dependence of the eigenvalues with e and s, we have that A(e, s) 
depends analytically on the variable e £ (0, oo) for any fixed s £ R^ and also on 
the variable s £ R^ for any fixed e £ [0, oo). In particular, A is differentiable 
in its second variable when the first is fixed. Moreover, A is continuous in both 
variables at the point (0, sq) for any sq £ R^ and A(0,0) = (0, 0). This is proved by 
proceeding exactly as in the proof of Proposition 12.11 cf. Section [3l to show that 


lim \Xk,e,s 

e \,0 


l^k,s 


= 0 . 


and 

(4.2) 


lim - 'yfc.slljLi(n) = 0 , 


for any fixed s £ R^ (in the proof one only has to change the reference metric go by 
the metric go,s)- Since this property does not depend on the particular way that 
the point (e, s) tends to (0,so)) we conclude that the map A is continuous on its 
domain. 

With the aim of using an implicit function theorem, we now need to evaluate the 
derivative of A with respect to s at the point (0, 0). To do this we first calculate the 
derivative of the Laplacian Ag, computed with the metric go,s, at the point s = 0, 
acting on the Neumann eigenfunction Vk- Noticing that AgUfe reads as 

^s'^k — 2 ^s: ^ V/e) /ifc/lg Vk 1 
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where V is computed with the reference metric go, a straightforward computation 
taking derivatives of this expression with respect to the parameter Si, i G {1,2}, 
and using the form of the conformal factor hg, shows that 

(4.3) Us, As(0)wfe = goiyipi, S/vk) + 2nki’iVk ■ 


Now we can take derivatives of the eigenvalue equation AgVk^s = with 

respect to Si (recall that the families of eigenfunctions Vk,s and fj,k,s were chosen to 
be analytic with respect to s), and evaluate at s = 0. If we multiply the resulting 
expression by Vk and integrate over the domain (with respect to the measure 
induced by go, which we omit for the sake of simplicity) we easily get 


(4.4) Ds.gk,s{,0) = - / VkDs.As{0)vk- VkA{Ds.Vk,siO ))-/ VkDs.Vk,s{,0) ■ 

J Q J Q Jq, 

To obtain this formula we have used that 

Noticing that the unit vector field i^s orthogonal to the boundary 912 with respect 
to the metric go,s is given by 

_ 1 

Us = hs V, 

with u the unit normal field of dfl computed with go, we infer that the Neumann 
eigenfunction Vk,s also satisfies the Neumann boundary condition with respect to 
the normal field v, that is, d^Vk^s = 0. In particular, taking derivatives with respect 
to Si in this expression we get d,^{Ds.Vk,s{0)) — 0. 

Accordingly, the llsiiifc,s(0) terms in Eq. (14.41) cancel after integration by parts, 
using that Avk = —gkVk- Finally, putting together Eqs. (Irm and (lia we obtain 
a Hadamard-type formula for the variation of the Neumann eigenvalue gk,s'- 


^Sigk,s{^) 



^'fe5o(VI/)i, Vufc) - 2g,k 



igki’i) . 




Here we have integrated by parts and used that d^tpi = 0, by assumption, in order 
to pass to the second equality. 

This equation and the fact that gi = g 2 = imply that the expression for the 
derivative ZIsA(0,0) is: 


n Ain nl = 1 r “ 4AiD2V’i) /n I'f (Ai /’2 - 4aid2?/’2) A 

* ’ 2\ J^vf{A^i - 4g,B2^i) v|(Ai/>2 - 4^d2i/’2) J 


In order to apply an implicit function theorem to the map A, in the follow¬ 
ing lemma we compute the rank of the matrix DgAiO, 0) for a suitable choice of 
the functions 1 pl,^|J 2 ■ The lemma is stated in terms of a function Vko that is any 
smooth extension to the whole manifold M of the Neumann eigenfunction Vkg of 
the domain H. 


Lemma 4.1. There exists a positive integer ko such that 

rank [L)sA(0 , 0)] = 2 
with the choice ipi =1 and ip 2 = Vko ■ 
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Proof. First notice that, with this choice, the functions tpi and '02 satisfy the Neu¬ 
mann condition = 0 on dPl. Since /q iifc = 1 by assumption (we are considering 
an orthonormal basis), the derivative I?sA( 0 , 0 ) reads as 


£>sA(0,0) 


If 4^d2 (^/co+ 4/iD2)/j^z;iz;fco \ 

2 V 4^d2 (m/co +4/iD2)/n^^2^fco j 


The problem then reduces to finding a Neumann eigenfunction such that 


[ {vf-V^)Vko ¥^0. 

Jn 

It is obvious that such eigenfunction exists because the set of Neumann eigenfunc¬ 
tions spans Lf{PL) and the quantity v\ — is not identically zero due to 

the fact that vi and V 2 are not proportional. This proves the lemma. □ 


We can now use the implicit function theorem m Theorem B] to obtain a 
constant eq > 0, an open neighborhood U of 0 GMf and a function R : [0, eo) —>■ U 
such that: 

• R{0) = 0, 

• A(e, i?(e)) = A(0, 0) = 0 for all e G [0, cq), and 

• R is continuous at e = 0, so in particular limex^o I?(e) = 0. 

We can define the desired conformal factor he as 

he ■= hR{e) , 

so that the continuity of R implies that limgx^o II ~ 1 ||c'(m) = 0 for any integer 1. 
With this definition of he we have that the metric and the eigenvalues Xk,e in 
the statement of the proposition are given by 

9e 9e,R{e) 5 '^fc,e,R(e) • 


Accordingly, since A(e, i?(e)) = 0 for e G [0, ep), we conclude that 

Al,e = A2,e = Ha2 

provided that e > 0 is small enough. The multiplicity of this eigenvalue is exactly 2 
because limex^p |A 3 ,e — ^al = 0 and /13 > 

Concerning the eigenfunctions Uk,e associated to the eigenvalues Xk,e, we observe 
that 

Uk,€ — 'f^k,e,R{e) 5 

and therefore Eqs. m and (14.21) imply 

1™ ll^i.e -'Cillffi(n) = lim Il?l 2 ,e-'i; 2 ||ffi(n) =0, 
thus completing the proof of the proposition. 


5. Higher dimensional analogs of the main theorem 

In Theorem o we have shown the existence of a Riemannian metric for any 
closed 3-manifold whose first nontrivial eigenvalue has multiplicity 2 (i.e., Ai = 
A 2 ¥ A 3 ) and a prescribed component of codimension 2 (that is, a knot) in the 
nodal set u“^(0) of the complex-valued eigenfunction u = ui-\- iu 2 . In this section 
we shall show that essentially the same argument enables us to construct a metric on 
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any closed d-dimensional manifold, whose first nontrivial eigenvalue has multiplicity 
m ^ 2 and a prescribed component of codimension m in the set 

However, a technical condition makes the statement considerably more involved 
in the general case. This condition is associated with the crucial requirement that 
the nodal set be structurally stable. In the situation covered by the main theorem, 
the structural stability follows from the important equation (12.211 . which plays a 
crucial role in the proof (see Remark 12.31) . The higher dimensional analog of that 
relation is the requirement that the intersection of the nodal sets uj"^(0 ),..., u“^(0) 
is transverse on S C that is, 

(5.1) rank^VMi(a;),..., Vum(T)^ = TO 

for all a; G S. 

We shall next state the higher dimensional analog of Theorem ll.il For this, let E 
be a compact embedded submanifold of a d-dimensional manifold M. Throughout 
we will assume that E has codimension to ^ 2 and is connected and boundaryless. 
Our result shows that E is a connected component of the transverse intersection 
of TO independent eigenfunctions corresponding to the same eigenvalue Ai of the 
Laplacian on M for some metric. For the transversality condition (lEID to hold, a 
topological obstruction on E is that its normal bundle must be trivial. Geometri¬ 
cally, this is equivalent to the assertion that a small tubular neighborhood of the 
submanifold E must be diffeomorphic to E x D"*, where D™ denotes the unit m- 
dimensional disk. Since a knot always has trivial normal bundle m, Theorem 11.11 
corresponds exactly to the case d = 3 and to = 2. The case to = 1, i.e. codimension 
one hypersurfaces embedded in M, was covered in [TO] . 

Theorem 5.1. Let T, be a codimension m eompaet submanifold of M with trivial 
normal bundle, to ^ 2. Then there exists a Riemannian metric g on M such that 
its first nontrivial eigenvalue has multiplicity to (i.e., Ai = ••• = Am )^m+i) 
and H is a connected component of the transverse intersection of the nodal sets 
uj"^(0)n- • •nM“^(0). Here {ui, • • • , Um} is any basis of the eigenspace corresponding 
to Ai. 

The proof of this result goes almost exactly as the proof of Theorem ll.il First, 
we construct a reference metric go on M such that 

9 o\q = dE + |dzp , 

where H is a tubular neighborhood of E, which can be identified with E x ID)'" 
because E has trivial normal bundle, and z £ are coordinates parametrizing 
the disk D™. We take the metric gg on E such that the first nontrivial eigenvalue 
of the Laplacian on E with g^, is bigger than the first Neumann eigenvalue /in™ of 
D™. Therefore, the Neumann Laplacian in fl with the metric go has eigenvalues 
Mo = 0, 

1^1 — * * * — /^m — 5 

and /ife > /iD”i for fc ^ to -I- 1. The eigenspace of /ii is then m-dimensional, with a 
particular orthogonal basis being 


Vk =r^ ™ {y/p,Dm r) Yi^k {6) = CmZk + 0{\z \‘^), 
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where Ji denotes the order Bessel function and Yi^k^ 1 ^ fc < m, are the first 
order spherical harmonics on The explicit expression of the constant Cm is 

not relevant for our purposes. From this formula for the eigenfunctions it is easy 
to check that the intersection 


is transverse and coincides with the submanifold S. The same holds for any other 
basis of the eigenspace corresponding to /ii. 

Step 2 is the same as in the proof of Theorem so that we can construct 
a 1-parameter family of smooth metrics gt = ftgo such that the corresponding 
eigenvalues Xk,t and eigenfunctions Vk,t behave when t — >■ oo as in Proposition 12. II 
In Step 3 we can prove a result analogous to Proposition [2^ The only difficulty is 
in choosing the functions tjji (see Section 0]) for the conformal factor 

U _ „2sitpi-\ - \-2Sm'lpm. 

Its — C , 


where s = (si, • • • , Sm) G R"*- The conditions that these functions must satisfy are 
that their normal derivatives on d^l must be zero, and the differential DsA{0, 0) 
has maximal rank, where the map A : R+ x R™ —>■ R™ is defined as: 


A(e,s) 


(^l,e,s , * * * , Xm,e,s ) if 6 0 , 

(Mi,s — Md™, • • • ) Mm,s — Mb™) ife = 0. 


The quantities Xk,e,s and fik,s are the eigenvalues corresponding to the metrics 
g^^s ■= hsgr, T = e“^, and go,s := ^s5o, respectively, with their associated or¬ 
thonormal eigenfunctions Ufc,£,s and Vk^s- By Kato-Rellich we can assume that all 
these quantities are analytic in the variables e > 0 and s G R™. 

The derivative DsA{0, 0) can be computed using a Hadamard-type formula as 
in Section m thus yielding: 

iDsA{0,0))ij = ^ vl{{d - 2)A^j - 


for 1 ^ i, j ^ m. To compute the rank of this matrix, we can prove an equivalent 
lemma to Lemma 14.11 by noting that the transversality of the nodal sets of Vk , 
1 ^ fc ^ m, implies that the functions are linearly independent and hence their 
span V is m-dimensional. This allows us to choose m linearly independent Neumann 
eigenfunctions {vk^, ■ ■ ■ ,Vk„} so that no function in their span is orthogonal to the 
subspace V (these eigenfunctions exist because the set of Neumann eigenfunctions 
is a basis of L^(n)). 

Therefore, if we set 

tpi = Vki 

for i G {1 ,... ,to}, where Vk denotes any smooth extension to the whole manifold 
M of the Neumann eigenfunction Vk, it follows that the vectors 


Wi := -- 


(^J ({d-2)g,ki+4:g,Bm'^ViVkm--- , J ((c^ - 2)Mfc„ + 4 md™) 


for 1 ^ i ^ m, are linearly independent, where each vector Wi is the z**' row of 
the matrix ZlsA(0,0). Indeed, if Wi were linearly dependent, there would exist a 
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nontrivial set of constants C K. such that 

m « 

/ v'lvk, =0 

i=i 

for all 1 ^ j ^ m. Hence, the function F := cipf would be a linear combi¬ 

nation of the functions vf that is orthogonal to each Vj-j ■ By the construction of 
the Neumann eigenfunctions Vkj , this would imply that F = 0 and from the linear 
independence of the set we would obtain that Ci = 0 for every 1 ^ ^ m, 

thus proving the claim. Notice that, as in the proof of Propositionwe can take 

V’l = 1- 

The rest of the proof goes exactly as in Section |4] using the implicit function 
theorem [H Theorem B] to prove a result analogous to Pror)osition l2.2l and Thom’s 
isotopy theorem in Step 4, so the details are omitted. 

Remark 5.2. Theorem o has the following striking corollary: since any exotic 
sphere of dimension m embeds smoothly in with trivial normal bundle m. 
we conclude that there is a metric on whose first nontrivial eigenvalue has 
multiplicity m and the exotic sphere is a component of the transverse intersection 
of the nodal sets uj"^(0) fl • • • fl u“^(0). Recall that an exotic sphere is a smooth 
manifold that is homeomorphic but not diffeomorphic to the standard sphere. 
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